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Abstract 

A reorganized perturbation expansion with a propagator of soft in- 
frared behavior is used to study the critical behavior of the mass 
gap. The condition of relativistic covariance fixes the form of the 
soft propagator. Finite approximants to the correlation critical ex- 
ponent can be obtained in every order of the modified, soft pertur- 
bation expansion. Alternatively, a convergent series of exponents 
in large orders of the soft perturbation expansion is provided by 
the renormalization group in all spatial dimensions, 1 < D < 3. 
The result of the e-expansion is recovered in the D — > 3 limit. 
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1. Introduction 



The critical behavior of various statistical systems near continuous phase tran- 
sitions has been the subject of intense investigations for several decades. Contin- 
uum quantum field theory turned out to be the most appropriate tool for these 
investigations. It connects critical phenomena with the fundamental features of 
renormalization and renormalization group. Nevertheless, the best numerical val- 
ues for critical parameters are obtained more often then not from lattice methods: 
high and low temperature expansions, Monte-Carlo investigations, etc. The rea- 
son for this dichotomy is obvious. Numerical methods on lattices are convergent 
procedures, which, as a function of parameters (or computer time), diverge only 
exactly at the critical point. In contrast to that, the critical behavior in continuum 
field theories is investigated through the e-expansion of Wilson and Fisher' 11 , which 
define only an asymptotic series for all values of the expansion parameter, e m 

In the current paper, we intend to introduce a new method for the investigation 
of critical behavior in quantum field theories. The method is based on a reorgani- 
zation of the perturbation expansion series, very similar to mass renormalization, 
with the intent of extracting key nonperturbative features from field theories. The 
purpose of the modification of the propagator is to soften the infrared behavior of 
the propagator, thus providing a finite approximant to the correlation critical ex- 
ponent in every order of the perturbation expansion. The condition of relativistic 
covariance fixes the form of the softened propagator up to a scale transformation. 
The form of the propagator is simply related to the correlation exponent, v. The 
series of approximants for u, obtained from successive orders of the soft pertur- 
bation expansion, is convergent, or at least has a convergent subseries. In this 
respect, our method is unique. 

The method will be demonstrated on the example of the 4 . In the next 
section, the modified perturbation expansion will be motivated and defined. Sec. 
3. will prove the existence of a soft propagator consistent with the requirement of 
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correct relativistic form of the energy of single particle states. The renormalization 
group approach will be applied to the theory defined by the modified perturbation 
expansion in Sec. 4. We will prove that an alternative series of exponents can be 
obtained in every spatial dimension, 1 < D < 3. Near D = 3 the leading order 
approximation reproduces the one loop result of e-expansion. In the final section, 
our results will be summarized. 



2. Modified Perturbation Expansion 

The stationary state Schrodinger equation for 4 field theory in D spatial 
dimensions is given by 

H \^) = lj d ° x (n 2 ^) + m2( t>( x ) + d^(x)d^(x) + = E\ij), 

(2.1) 

The investigation of a field theory in Schrodinger representation is easiest using 
the representation 

n(x) = (2.2) 

i d(p(x) 

for the canonical momentum operator. Then one can choose to be a functional 

f3l 

of the classical field 0. This corresponds to the coordinate representation of 
quantum mechanics. 

We will investigate the equations for the ground and first excited (single parti- 
cle) states simultaneously. Our aim is to extract information concerning the critical 
behavior of the renormalized mass. We will denote the energy eigenvalues of the 
ground state and the first excited state by and E^> = E(q) + E^°\ respec- 
tively. The first excited state of the system is a single particle state. The excitation 
energy of a single particle state must have the correct relativistic form 

£ 2 (q) = q 2 + m 2 , (2.3) 

where the notation m is used for the physical mass. 
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As usual, we rearrange the Hamiltonian by subtracting a mass renormalization 
term, 



from the kinetic part and adding it to the interaction part. Then, with an ap- 



D < 3. We define the subtraction scheme more precisely below. 

The Hamiltonian perturbation expansion (or old fashioned perturbation ex- 
pansion) is completely equivalent to the covariant perturbation expansion. After 
collecting terms in every order of the perturbation expansion, one can introduce 
the extra po momentum components and get an expression for the squared energy 
function E 2 (p) as 



where £(po, p) is the sum of self-energy corrections. The self-energy correction is 
covariant, thus, it depends onp 2 = pjj + p 2 om y> but for future reference we indicate 
its dependence on the energy component independently. Furthermore, if m is the 
renormalized mass, the self-energy correction must be subtracted as S(m, 0) = 0. 
Then, of course, every term of the covariant perturbation expansion of S, being 
momentum independent, vanishes for all p. The reader may ask at this point, why 
do we go through this elaborate procedure to define an identically zero quantity? 
The answer is that we will use a similar procedure later, in which the vanishing of 
E will not be trivial at all, and at that point we will draw on the analogy to the 
trivial excercise above. 

The information concerning the critical behavior extracted from the perturba- 
tion expansion is of a divergent nature. The root of this problem is shown clearly 
by the m dependence of the terms of the perturbation expansion of 5m 2 . Dimen- 
sional considerations show that in nth order of the perturbation expansion, after 




(2.4) 



propriately chosen 5m 2 = m\ — m 2 , the perturbation expansion becomes finite for 



E 2 {p) = p 2 + m 2 + £(R),p) 




Po= p' 2 +m 2 
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the removal of a term diverging when the cutoff, A — > oo, but independent of m, 
one is left with a term depending on the renormalized mass as 



5m 2 



c n g /I 



,n n{3-D) 2-(3-D)n 



(2.6) 



where c n is a dimensionless constant. For every 1 < D < 3 and for sufficiently 
large n, the power of m becomes negative. On the other hand, the dependence of 
the 'exact' 5m 2 on m is related simply to the critical behavior 



where A > \ jv and v is the correlation critical exponent and m 2 (A) is the critical 
value of the squared bare mass. Thus, the m — > behavior of (2.6) has nothing to 
do with the real critical behavior. 

Our thesis is that one does not obtain the correct critical behavior because some 
crucial non-perturbative features of the theory are overlooked by a straightforward 
perturbative expansion. Our aim is to incorporate some non-perturbative features 
into the expansion such that the critical behavior of the renormalized mass can 
be extracted. The procedure we intend to employ is very similar to the additive 
renormalization of mass. 

There are examples of dramatic improvement of the convergence of pertur- 
bation series by reorganizing the perturbation series. Halliday and Suranyi were 
able to obtain convergent perturbation expansions for the anharmonic oscillator in 
quantum mechanics by an appropriate reorgnization of the series! 4 " 01 Following the 
spirit of these works, we will alter the propagator by adding a term 



to the kinetic part of the Hamiltonian and subtracting it from the potential part. 
By this rearrangement, we are not able to make the perturbation series convergent 



5m 2 = 



m 2 c (K) + cm* + 0(m x ), 



(2.7) 




(2.8) 
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for non- universal quantities like the coefficient of the term m 1 /^, but we are able 
to extract convergent approximations for the exponent v. 

At this point, the form of the corrected propagator is undetermined. Using an 
arbitrary propagator, however, changes the perturbation series in an uncontrolled 
manner, thus, the series expansion for E(p) will not have the correct relativistic 
form (2.6). We will show in the next section that requiring the correct relativistic 
form for L 2 (p) in every order of the perturbation expansion fixes the form of A(p), 
up to an irrelevant overall multiplier. 

It will be convenient to use the D + 1 dimensional, 'relativistic' form of the 
perturbation expansion. Then the propagator will take the form 

D(P,P0) = ~ 2 1 r— - , (2.9) 

Po~u (P) 

where 

u 2 (p) =p 2 + m 2 + A(p), (2.10) 
while the interaction term will be 



(2.11) 



£int = \f d D+1 x [fy\x) +<ymV(:r)] - \ J ^ +1 pA(p)|0(p)| 2 
= L[ nt + Lg m 2 + La, 

where for convenience we have given La in momentum representation. 

It is easy to see, starting from the D dimensional, Hamiltonian form of the 
theory, that in the language of the D + 1 dimensional 'relativistic' theory, the 
quantity we must calculate is 



L 2 (q)=cu 2 (q) + S(g ,q) 



(2.12) 

^=w 2 (q) 



In contrast to the similar expression obtained in the standard perturbation expan- 
sion, (2.5), one does not trivially obtain the correct relativistic form for the single 
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particle energy. Requiring the correct relativistic form (2.6) implies the nontrivial 
integral equation 



A(q) + E(g ,q) 



= (2.13) 

gg=q 2 +m 2 +A(q) 



for A(q). 1 



There is a trivial solution of (2.13), A(q) = 0, reducing (2.13) to (2.5). We 
will prove in the next section that a nontrivial solution for A(q) exists as well such 
that in every order of the modified perturbaton expansion (2.13) is satisfied and 
the correct dispersion relation for the energy is maintained. 

3. Existence of the Modified Expansion 

The perturbative terms become bounded functions of the physical mass, m, 
unlike (2.6), only if the infrared behavior of the propagator is softened. Thus A(p) 
should have a leading order infrared behavior 

A(p) ~ |p| 2a e 2 " 2Q , (3.1) 

where < a < 1 and £ is a scale parameter of the dimension of a mass. Then the 
p dependence of the propagator is dominated by A(p) in the interval J, 

I:£»\p\»m. (3.2) 

The euclidean propagator of the 'relativistic' perturbation theory has the form 



1 In a previous work, wc obtained results for critical exponents, identical to the ones ob- 
tained in this paper, bj 
interpretation in Sec. 5. 



tained in this paper, by a different interpretation of (2.12)! 61 We will comment on that 
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In interval /, the propagator can be approximated by 



1 



(3.4) 



Po + IP 



where we absorbed the scale factor into momentum p (in other words, we set 



For the time being, we will omit term La of (2.11) in the perturbation expan- 
sion. We will prove later that its inclusion would not change our results for the 
critical exponents. 

It will be sufficient to study (2.13) in interval / to fix the form of A(p). We 
can argue that if we wish to extract the leading behavior of the terms of our 
expansion while the external momentum is in interval J, then it is self consistent 
to use propagator (3.4), instead of the complete form (3.3). That is certainly so if 
the loop integrals are dominated by contributions in which all internal momenta 
are in interval / as well. 

An n loop contribution generated by such propagators has the behavior |q| a ", 
where a n = 2a + n(D — 3a). This can be seen if we realize that in a 4 theory there 
are 2n — 1 propagators, each contributing |p| _2c \ while the momentum integrations 
over p andpo contribute \p\ Dn and |p| nQ , respectively 1 . Using similar arguments, it 
is easy to see that n loop vertex corrections have the overall momentum dependence 
|q| 6n , where b n = n(D — 3a). Now a necessary condition for the finiteness of self 
energy corrections is the finiteness of vertex corections. Those are finite only if 



D 

a > «o — — • 



(3.5) 



We will restrict ourselves to such values of a in the future. 



1 Note that po — 0(|p| a ) in loop momenta. 
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Suppose now that in a certain loop integral the loop momentum is in the region 
|p| > A 4 ) but the momenta external to the loop are in /. Then the loop is effectively 
shrunk to a point vertex with at least four external legs (self energy insertions will 
be discussed below). If the loop contains k propagators (k > 2), then the overall 
behavior of the diagram has been increased by a power of 2ka — a — D > 0, in 
view of the constraint (3.5) on a. Then we do not obtain a leading contribution, 
since |q| << £ = 1. 

If the external momentum of a self energy insertion is in /, but internal mo- 
menta are not, then the correction is even smaller. This can be seen by observing 
that self energy corrections are subtracted at pjj — P 2 — ( m2 can be neglected 
in /) and the largest contribution must be proportional to Pq. Then observing 
that Pq = 0(|p| 2a ) >> |p| flfe , for all k > 1, we can see that again we obtain a 
non-leading contribution. Thus we have proved that if the external momentum is 
in interval J, then propagator (3.4) should be used. 

In iVth order of the modified expansion, the infrared behavior of the diagrams 
is controlled by power a n = 2a + n(D — 3a). If a n > 0, diagrams are infrared 
convergent, even at vanishing external momentum. We can easily see that the 
subtraction at qo = q = makes the contribution finite as long as b n = n(D— 3a) < 
0, where ultraviolet divergences will appear. Then the subtracted n loop integral 
is finite in interval J n 

D D 

— 7T > a > — 
3- ^ 3 



Jn ■ 2 > a> —■ ( 3 - 6 ) 



n 



In the subtraction term, we set go = and q = 0. The subtraction results in 
differences having large momentum behavior 



|2a |n|2 a lnl 2a— 2 

Ml I I 1 I 



alql 2 - 2ap -q+2a(a - 1) ^'^ + ... 

P 



, (3.7) 



where p is a loop momentum. Since the leading correction term, — 2ap ■ q, cancels 
after angular integrations over p, the ultraviolet power behavior is reduced by two. 
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The subtraction in an expression of the form (po ~ Qo) 2 ~ Po = Qo ~ ^Polo reduces 
the power behavior by ~ |p| 2a only 1 . In other words, the ultraviolet behavior of 
the subtracted self-energy correction term is controlled by the power a n — 2a = b n . 
This is, of couurse, the same power as the one controlling the ultraviolet behavior 
of the vertex part. Our assertion concerning the convergence of the term of the 
expansion in interval J n , (3.6), has been proved. 

3.1. Existence of a nontrivial solution 

Next we show that a value for a G J n can be found such that the leading 
term in (2.13) vanishes. Non-leading order terms that might appear in (2.13) and 
non-leading terms of A(q) will be considered later. Suppose we work in n loop 
order of our approximation scheme. Then if |q| is in interval J, the largest order 
(n loop) term dominates. This is so, because the series of exponents a& satisfies 
a k > a k+l- Since |q| 2Q << |q| a ™, solution (3.1) of (2.13) is possible only if the 
coefficient of the term |q| an vanishes. 

First, in a given order of the loop expansion all euclidean diagrams are real and 
have the same sign. The subtraction (mass renormalization) term of every diagram 
becomes infrared singular at a n =0. Since the subtraction term is always negative, 
the n loop term tends to — oo if a — > D/(3 — 2/n). On the other hand, it is also 
easy to see that the n-loop contribution tends to +oo at a — > D/3, the point of 
ultraviolet divergence. The contributions are analytic functions of the variable qfa. 
In other words, we can write a dispersion relation for the subtracted contribution 

oo 

J L % + z z 

where the positive definite discontinuity has an asymptotic behavior f(z, q) ~ 
z a n /2a_ jsjote that the constraint < a n < 2a insures the convergence of (3.8). 

1 Note that |p| 2Q >> p 2 , if |p| e /. 



(3.8) 
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Dimensional considerations show that |q| can be scaled out from (3.8) to give 

oo 

\q\ aN J dz 

Furthermore, according to our previous discussions, the leading asymptotic behav- 
ior of the n loop term comes from the term in the denominator of the integrand. 
Thus, at large values of z, we can also set the momentum 1 equal to zero in the 
argument of the spectral density. Finally, according to our prescription, we have 
to set g 2 = — q 2 — m 2 — A(q) ~ — |q| 2a , for the Euclidean go- Then we obtain 
the following expression for the large z contribution, responsible for the ultraviolet 
divergence 

oo 

|qr| dzf(z,0)—^—-. (3.10) 

(3.10) shows that the integral indeed tends to +oo when a — > D/3. We can 
conclude that the n loop contribution to the subtracted self energy diagram has 
a zero at some a n G J n . At a = a n , the coefficient of the leading power, |q| a ™, 
where a n = 2a n + n(D — 3a n ), vanishes. 

We turn now to the discussion of non-leading order terms in A(q). Though by 
setting a = a n , the n-loop contribution vanishes, the correct dispersion relation 

r\ k 

for E 2 (q) requires that the coefficients of the powers |q| an (A;-loop term), where 
a\ = 2a n + k(D — 3a n ), k < n should also vanish. For k > 0, these powers are 
also smaller than 2a and give contributions much larger then |q| 2a in the region 
|q| << 1, so they should also be cancelled in (2.13). It is easy to see, however, that 
a non- leading term of A(p) in the propagators of the loop diagram can cancel these 
contributions as well as the term A(q) appearing in (2.13). Thus A(p) should have 
the form 

A(p) = ^^Q|p| 2a+ ^ M ). (3.11) 

1=0 

Then an expansion around the leading contribution co|p| 2a leads to correction 
terms to the n loop contribution, which have exactly the same behavior as the 



q 2 /\q\ 2a + z z 



(3.9) 
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n — 1, n — 2, etc. loop terms. An appropriate choice of the coefficients q in 
(3.11) cancels these contributions. It is amusing to notice that in large orders 
of the perturbation expansion the coefficients q decrease very fast with n, e.g. 
ci ~ cjj/n 2 . 

Finally, let us return to the neglected term, of the perturbation expansion. 
It is easy to show that the contributions of L/± are similar to those of fc-loop terms, 
k < n, of the perturbation expansion. Indeed, if we work in n loop order, an 
insertion A(p) on a propagator would substitute a four point vertex, lowering 
the power of g by one, and would substitute one of the propagators D(qo,p) by 
D(qo, p)A(p)D((7o, p) = 0(D(pQ,p)). In other words, we obtain a contribution 
with the same power of g and the same low momentum behavior as the n — 1 loop 
term. We know, however, that the n loop term provides the dominant power of 
the expansion, thus the equation determining a is not affected by the term 
of L', defined in (2.11). These non-leading terms can also be cancelled in (2.13) 
successively, by non-leading terms of the asymptotic expansion of A(p). 

3.2. Critical exponent 

The correlation critical exponent v relates the behavior of the unrenormalized 
mass and the renormalized mass 

ml = m 2 + 5m 2 = m 2 + /(A) + cm* , (3.12) 

where /(A) is the cutoff dependent critical point for the square of the bare mass 
and c is a constant. The exponent a n has a direct connection with the nth order 
approximant of v. One is able to extract v if one investigates the behavior of 
the terms of the modified expansion for 5m 2 in the transition region, |q| ~ m. 
In that region, the mass term cannot be neglected and we have to use the form 
Pq + \p\ 2a + m 2 for the propagator. 
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Due to the infrared convergence of the integrals in interval (3.6) we can separate 

7 k 

terms of the form c^g h OLn from the mass renormalization constant, such that the 
remainder is a finite function m. Then we can scale out the mass, m, from the terms 
of the perturbation expansion with the simultaneous substitutions p — > m~p and 
po — > mpQ. Using (3.12), we obtain the following expression for v n , the estimate 
for the critical exponent v in n-loop approximation, 

— = 2-n(3 ). (3.13) 

Equation (3.13) shows that the dominant contribution (smallest power of m) 
comes from the leading order, n-loop term, since for k < n we have 2 — k(3 — 
D/a n ) > 2 — n(3 — D/a n ), provided a G J n of (3.6). Furthermore, it also follows 
from (3.6) that the value of v is between \ and oo, as it indeed should be. 

The natural method for finding approximants to v would be to find values for 
a n at as high n as possible. Indeed, in a previous paper we found the critical 
exponent for D = 1 and D = 2 in the two-loop approximation, i/ 1 ) = 1.013, 
z/( 2 ) = 0.713, in reasonable agreement with known values. 161 

There is a puzzling question concerning our approximation scheme. Since a n 
is obtained from n-loop diagrams only, completely different from the (n — l)-loop 
diagrams used for obtaining a n -i, h is hard to see the relation between subsequent 
approximants to v, v n and v n -l- It could happen that function of n, 

oscillates between \ and oo, instead of being convergent. A bounded infinite series 
always has at least one convergent subseries, but it may have several that converge 
to different limits. As is shown by (3.13), a smooth limit to the correct value of 
v is possible only if we can prove that in large orders of the loop expansion the 
equation determining a becomes a function of p = n(3 — D/a) only, thus making 
the existence of a nontrivial limit p — > p c for n — > oo possible. Then, of course, 

2— p c 
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The simplest method for the investigation of large orders of the perturbation 
expansion is the renormalization group. Observe that the modified loop expansion 
is equivalent to a perturbation expansion in an euclidean field theory with action 1 



Similar field theories have been investigated in the past. The closest example 
is that of Fisher and Grinstein who found that certain superconducting systems 
may be described by Hamiltonians similar to (3.14). They, however use a 'double e- 
expansion' in the variables e = 3 — D and e' — 1 — a. Our intention is to investigate 
field theory (3.14), at fixed D near the critical point a = D/3. In other words, we 
wish to perform an expansion in 5 = 3a — D. 

We will show in the next section that the Callan-Symanzik equation may be 
used to find the coefficients of terms, singular in 5, in large order contributions to 
Green's function. This is quite analogous to the calculation of leading logarithms 
in perturbation theory. We will also find that up to a multiplier, independent of the 
momentum, the large order contribution to the self energy correction becomes a 
function of the variable p = 3n5/D. Higher and higher orders of the loop expansion 
provide a power series of higher and higher order in variable p. Finding the zero 
of the self energy correction as a function of p provides an alternative method for 
the determination of v. One outstanding feature of the power series expansion in 
p is that in contrast to the e-expansion, it has a finite radius of convergence. 



1 For the sake of simplicity, we mix coordinate space and momentum space representations 
for </>. The action could be written in either of these representations without any difficulty 




dxod D x<f)(x) 4: + I dpQd D p(/)(p)(f)(~p), 



(3.14) 
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4. Renormalization Group 



4.1. Callan-Symanzik equations 

The renormalization of the field theory, defined by action (3.14) requires wave 
function and coupling constant renormalizations, besides the additive mass renor- 
malization we discussed earlier. The renormalized and bare amputed Green's func- 
tions have the following relation with each other 

r (4) (pi,P2,P3,P4,#,/i) = Z- 2 (g,5)T i Q\p 1 ,p 2 ,P3,P4,go,S), (4.1) 

and 

T^(p,g^) = Z- 1 (g,5)T^\p,g ,5), (4.2) 

where 5 = 3a — D and \i is the renormalization scale, defined through the relation 
of the renormalized and the unrenormalized charges 

g Q = 9l i s Z g {g,S). (4.3) 

Pi denote both the spatial, D dimensional part of the momenta, and their timelike 
part as well. is the inverse propagator, the main subject of our investigation. 

The derivatives of the left hand sides of (4.1) and (4.3) with respect to 1/5 
should vanish as 5 2 , when S — > 0, because their first nonvanishing 5- dependent 
terms are proportional to 5. Thus the unrenormalized Green's functions satisfy 
the following equations 

r^(p h ...p 4 ,g,5)=0(5 2 ), (4.4) 
and 

T^\p,g,5) = 0(5 2 ), (4.5) 

where we dropped the subscript of the coupling constant. Strictly speaking, it is 
not exactly r^ 2 ) that satisfies (4.5), as it will soon become obvious. The relation 
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between the standard cutoff dependent form of the Callan-Symanzik equations and 
the above forms is given by the two possible limits 

ln(A/ „ _ hzS^l _ i. (4 . 6) 

While the cutoff scheme is obtained if one takes the limit 5^0 first, the di- 
mensional regularization scheme is obtained if the limit A — > oo is taken first. 
The beta and gamma functions can be determined from (4.4) and (4.5), taken at 
specific momentum values. In the 5^0 limit they coincide with the standard 
definitions. 

The renormalization group functions have the following perturbation expansion 

P(g) = fog 2 + fog 3 + - (4.7) 

and 

70?) =72£ 2 + 733 3 + - (4-8) 



4.2. Calculation of self energy corrections in large orders of 
the soft perturbation expansion 

The two point amplitude can be written in the form 

if \ P0 , P , g, 5)=pl + |p| 2a + p§S(po, P, g, 5), (4.9) 

where we explicitly indicated the separate dependence on po and p. The important 
property of T^ 2 ) is that the coefficient of |p| 2a does not get renormalized. When 
the subtraction at zero momentum is performed in the loop integrals, the leading 
terms in the difference are of the form (po — go) 2 + (po + <?o) 2 — ^Qq = Po and 
(P - q) 2a + (P + q) 2a - 2q 2a ~ p 2 << Pq, if p << £■ That of course shows that 
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indeed, the coefficient of p 2a does not get renormalized, and we should drop the 
second term from the right hand side of (4.9), when we substitute into (4.2) or 
(4.5). 

We expand the function S in g as follows 

S(po, P, 9,5) = Y1 {jZ ^T L Qn{^ P), (4.10) 

n 

where the function Q has the series expansion 

71-1 

The leading coefficient q n (p) = 1, by the normalization we chose for Q n , while 



7 n— 1 

m 

subtraction scheme 



q™ 1 (p) is the finite part of the n loop amplitude, divided by (— fa) 71 , in the minimal 



Callan-Symanzik equation (4.5) and (4.10) imply the following equation for Q n 

= (n-l)[Q n -Q nr . 1 ]-6^ 
+ §(n - 2)5Q n - 2 - §(n - 3)Q n _ 3 + ... ( 4 - 12 ) 

Now Q n is a slowly varying function of n. Then, in leading order of n and fixed 
x = 5n one obtains the following differential equaton for Q n 

n dQ n dQ n 03 n 

= n-^-5- Br + -^xQ n . (4.13) 
Solving (4.13) one obtains the solution of (4.12) as 

Q n = F{x) exp j^xlognj [1 + 0(l/n)], (4.14) 
where F(x) is to be determined by the inital conditions. 
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Note that all the momentum dependence is concentrated in F(x). Nevertheless, 
since F(x) is dimensionless, it becomes independent of the momentum components 
after the substitution po = |p| Q . The zero of F(x), as a function of x determines 
the value of a n . As x = Dp/3, it provides solutions at fixed p, as required by the 
existence of a nontrivial asymptotic value for the critical exponent. 

The initial conditions fix the form of F(x). In fact iterating the expanded form 
(4.12), with @3 = 04 = ... = 72 = 73 = ... = we obtain 



As we mentioned before, the coefficient q^ + i is the finite part of the m loop term, 
divided by (— P2) 171 , in the minimal subtraction scheme. 

It is well known that in a 4 field theory the asymptotic m dependence of m 
loop contributions to Green's functions T n ~ (—l) m m\a~ m m k b. Consequently, 
expansion (4.15) of F(x) has a finite radius of convergence, a. The improved 
convergence of series (4.15) compared to that of the e-expansion is due to the 
extracted multipliers n m , contained in the power x m . 

The finite radius of convergence, and the fact that the singularities of the 
expansion series (4.15) are at x < (the series is alternating) allows a conformal 
transformation of the variable x such that all values < p < 2 lie in the domain 
of convergence. 

Unfortunately, due to the complicated form of the propagator of the effective 
theory we have not been able to determine parameter a, controlling the radius of 
convergence yet. We know, however that F(x) must have a zero between < p < 2. 




n-l 



(4.15) 
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4.3. Critical exponents in first approximation 



The leading order term of integral equation (2.13) vanishes if a is the root 
of the subtracted n-loop contribution. As we have seen it earlier, at large n this 
is equivalent to finding the root of the n loop contribution as a function of x = 
n(3a — D) = n5. Then v = 1/(2 — 3x/D). Since the n-loop contribution is 
proportional to F(x), we have to find a root of F(x). Then a series of critical 
exponents, Uk, alternative to the series obtained in Sec. 3, can be obtained by 
approximating F(x) by its terms up to x k . 

In particular, for k = 1 we obtain x\ = —l/q^, where q\ is the ratio of the 
minimally subtracted finite term and of the pole term of the subtracted two loop 
self-energy diagram. Let us determine q\. It is simpler to use old- fashioned per- 
turbation theory (the amplitude integrated over timelike momentum components). 
The subtracted two loop amplitude, up to irrelevant multipliers, is 



(2)= [cPpid^i 1 

J Pf P2 l[p? + Pf + (i-Pi-p 2 )1 2 -i 



[P? + Pf + (Pi + P2) 



al2 



(4.16) 



In (4.16), as required, we substituted = — q 2a for the components of the external 
momentum. Furthermore, we factored q 2I) ~ 4a out. i is a unit vector. 

It is comparatively easy to separate the pole term from (4.16), in variable S. 
We obtain 



r (2) = 
pole 



5 



d£ / dQ D 



£2D/3- 



[1 + e C/3 + (1 + £ 2 + 2£ cos 9) D l & f ' 



(4.17) 
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The finite part of (4.16), taken at a = D/3 is 



r (2) f d D Pl d D P2 \ 1 

1 finite / D/3 D/3 | r D/3 . D/3 . ( z \D/3l2 i 

J Pi p 2 Mpi + P2 + (i - pi - V2) D,i Y ~ 1 

1 ^(|P1|-1M|P2| 



[pf /3 + P ? /3 + (pi + p 2 )^ 3 ] 2 [pf /3 + P ? /3 + (pi + p 2 ) / 3 ] 4 

+ ~T / ^ / dQ ^ 2D/3 ~ 1 { [1+^/3 + {1 + e + 2 ^ cose) D/6 ] 4 


^D/3, \ nx + 2 ln(l + x 2 + 2x cos 0) (1 + x 2 + 2a; cos 9) D / 6 
[1 + e c/3 + (1 + £ 2 + 2£ cos 6) D l & f ' 

(4.18) 

(2) (2) 

We have evaluated the ratio x = — ^pole/^&Ate ^ or D — 2. We obtained 
x = .402 , resulting in the prediction for the critical exponent, v = .717, very 
similar to the one obtained by our first method. Furthermore we examined the 

(2) 

limit D — > 3 as well. The first term of Tg n ; ite develops a pole at D = 3. Using the 

(2) (2) 

standard notation e = 3-D, one obtains rp ole /rg n ; ite = — e/3, and v — 1/(2 — e/3), 
in agreement with the leading order result from e- expansion . 



5. Summary 



A modified perturbation expansion was introduced in this paper. The expan- 
sion has the virtue of providing a potentially convergent series of approximants for 
the critical exponent v. The subsequent approximations are determined by subse- 
quent orders of a modified perturbation expansion, in which the infrared behavior 
of the propagator is softened. The condition of relativistic invariance fixes the 
propagator uniquely. 

The results of the first part of this paper are identical to those of a previous 
work of ours) 61 in which we arrived at the approximants to the critical exponent 
taking a different route. In that paper we obtained a slightly different form of 
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(2.12), 



E 2 (q) = q 2 + m 2 + S(g ,q) 



(5.1) 

q 2 =E^q) 



which was interpreted as an integral equation for E 2 (q). Then a solution of non- 
relativistic infrared behavior, E 2 (q) ~ cq Q , was obtained, where a coincides with 
the value of a obtained in the current paper. In Ref. [6] it was assumed that the 
self-consistency of (5.1) requires the vanishing of coefficient c in infinite order. 

The results of the second half of our paper are even more interesting. The mod- 
ified perturbation theory defines an effective field theory in the infrared domain. 
Large orders of that effective field theory can be investigated through the renor- 
malization group. The k-loop approximation to the solution of the renormalization 
group equation was shown to provide a series of finite radius of convergence, 
the zero of which also provides an approximation to the correlation exponent. 

The prescription, we provided for the determination of the correlation exponent 
in n loop order uses the n loop terms alone. This is the very circumstance which 
makes the universality of the critical exponent (its independence from the coupling 
constant) possible. In the e-expansion approach the dependence on the coupling 
constant is eliminated by the limit g — > g c (e), where g c (e) is the fixed point of the 
renormalization group transformation. Here we approach the critical surface at 
an arbitrary value of the coupling constant, by adjusting the bare mass. Thus, v 
can be independent of g only if in every order it is determined by the highest loop 
contribution alone. 

Finally, we discuss the extension of the results of this paper to the O(N) 
model. The only difference a calculation in the O(N) model would make is that 
the contribution of individual diagrams would be multiplied by an N dependent 
multiplier. There is only one diagram contributing to the equation determining 
a n in two and three loop order (Note that we have chosen the subtraction scheme 
in which momentum independent insertions on propagators cancel). Thus only in 
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the four or higher loop orders (where there are four or more topologically different 
diagrams) would we get a iV-dependent v. This, is certainly not a fatal flow, 
since the dependence of v on iV is weak. As far as the N dependence of the 
critical exponent is concerned the situation is similar in the renormalization group 
approach. This is a little more surprising, because for N = 1 (but not for other N) 
at D ~ 3 the result of the e-expansion is reproduced in the two loop approximation. 
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